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ABSTRACT 

We investigate the origin of the shape of the globular cluster (GC) system scaling pa- 
rameters as a function of galaxy mass, i.e. specific frequency (Sjv)i specific luminosity 
(Sx)> specific mass (Sm), and specific number (T) of GCs. In the low-mass galaxy 
regime (My >— 16 mag) our analysis is based on EST /ACS observations of GC popu- 
lations of faint, mainly late- type dwarf galaxies in low-density environments. In order 
to sample the entire range in galaxy mass (My = — 11 to —23 mag = 10 6 — IO^Lq), en- 
vironment, and morphology we augment our sample with data of spiral and elliptical 
galaxies from the literature, in which old GCs are reliably detected. This large dataset 
confirms (irrespective of galaxy type) the increase of the specific frequencies of GCs 
above and below a galaxy magnitude of My ~ —20 mag. Over the full mass range, the 
Sl— value of early- type galaxies is, on average, twice that of late-types. To investigate 
the observed trends we derive theoretical predictions of GC system scaling parameters 
as a function of host galaxy mass based on the models of Dekel & Birnboim (2006) in 
which star-formation processes are regulated by stellar and supernova feedback below 
a stellar mass of 3x 10 10 A1©, and by virial shocks above it. We find that the analyt- 
ical model describes remarkably well the shape of the GC system scaling parameter 
distributions with a universal specific GC formation efficiency, n, which relates the 
total mass in GCs to the total galaxy halo mass. Early-type and late-type galaxies 
show a similar mean value of 77 = 5.5 x 10~ 5 , with an increasing scatter towards lower 
galaxy masses. This can be due to the enhanced stochastic nature of the star and 
star-cluster formation processes for such systems. Some massive galaxies have excess 
n values compared to what is expected from the mean model prediction for galaxies 
more luminous than My ~ —20 mag (Ly > 10 10 L Q ). This may be attributed to a 
very efficient early GC formation, less efficient production of field stars or accretion of 
predominantly low-mass/luminosity high— 77 galaxies, or a mixture of all these effects. 
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1 INTRODUCTION 

One of the very first stellar systems to form in the early 
Universe are old globular clusters (GCs), which are observed 
in galaxies of all morphological types. Globular clusters are 
massive agglomerations of gravitationally bound stars, the 
majority of which formed almost simultaneously out of gas 
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with similar chemical composition. As fossil records of the 
first star formation episodes of their host galaxies, the distri- 
bution of their integrated properties (age, mass, metallicity, 
structural parameters) as well as the general properties of 
all GCs in a galaxy (total numbers, spatial and dynami- 
cal distributions) hold important clues to the initial physi- 
cal conditions at which they have formed and evolved. For 
that reason, the global properties of globular cluster sys- 
tems (GCSs) have long been recognized as promising tools 
to study the major galaxy star formation episodes and to 
serve as observational constraints to differentiate between 
various models of galaxy formation (Kissler-Patig 2000; van 
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den Bergh 2000; Harris 1991, 2003; Brodie & Strader 2006, 
and references therein). 

One of the most commonly used parameters to describe 
GCSs is the specific frequency (Sn), i.e. the number of GCs 
per unit galaxy luminosity. In essence, Sn measures the for- 
mation efficiency of GCs relative to field stars. GCs and 
field stars are linked through the dissolution of star clusters 
due to various mechanisms (e.g. tidal shocks, cluster relax- 
ation), which shapes the initial power-law to the observed, 
present-day Gaussian globular cluster luminosity function 
(GCLF, see e.g. Fall & Zhang 2001; Goudfrooij et al. 2004, 
2007; Gieles et al. 2006; McLaughlin & Fall 2008; Kruijssen 
& Portegies Zwart 2009; Elmegreen 2010, and references 
therein). The Sn parameter was initially introduced by Har- 
ris & van den Bergh (1981) as a measure of the richness of 
GCSs in elliptical galaxies. Since then Sn has been applied 
by numerous studies to galaxies of different morphological 
types from early-type, quiescent ellipticals to actively star- 
forming, late- type spirals, irregulars, and interacting/merger 
galaxies, covering the entire range in galaxy mass (from gi- 
ants to dwarfs) and environments (from galaxies in dense 
clusters to such in loose groups and in the field). Wide-field 
ground-based and deep Hubble Space Telescope (HST) stud- 
ies show that the Sn value varies greatly among galaxies, 
particularly among the most luminous ellipticals and low- 
mass dwarf galaxies (e.g. Harris 1991, 2001; Wehner et al. 
2008; Peng et al. 2008). Spiral galaxies tend to show much 
less scatter in their Sn, with values in the range 0.5 — 2 (Ash- 
man & Zepf 1998; Goudfrooij et al. 2003; Chandar et al. 
2004; Rhode et al. 2007). The general trend observed for 
early-type dwarf elliptical (dE) galaxies is that their aver- 
age Sn value increases with decreasing galaxy luminosity 
from a few to a few tens (Durrell et al. 1996; Miller et al. 
1998; Miller & Lotz 2007; Peng et al. 2008). It was shown 
recently that a similar behavior of Sn holds for late-type 
dwarfs (dlrrs, Seth et al. 2004; Olsen et al. 2004; Sharina 
et al. 2005; Georgiev et al. 2008; Puzia & Sharina 2008). 
Conversely, in the high-mass galaxy regime, a significant in- 
crease of the Sn values is observed for the most luminous 
galaxies (gEs and cDs, e.g. Rhode et al. 2005; Peng et al. 
2008; Harris et al. 2009). As the stellar populations (colours 
and integrated light spectral indices) of the latter galaxies 
are not significantly different from galaxies 1-2 mag fainter, 
such an upturn in the Sn distribution implies that the as- 
sembly history of the most massive galaxies must have been 
different from those of less massive galaxies. 

The overall trend for Sn is that with increasing galaxy 
luminosity from My ^ —11 mag to Mv ^ —20 mag, Sn 
decreases from the range ~ — 100 to Sn ~ 0.5 — 3, respec- 
tively. For galaxies more luminous than Mv — —20 mag the 
Sn increases again (to ~ 10). 

To explain the increasing fraction of GCs over field stars 
in the most massive galaxies several studies suggested that 
the mass in GCs is proportional to the total gas mass sup- 
ply (West et al. 1995; Blakeslee et al. 1997; Blakeslee 1999). 
In particular, McLaughlin (1999) investigated the Sn be- 
havior for the entire early-type galaxy mass range, from gi- 
ants to dwarfs. He showed that the enhanced Sn for the 
most massive galaxies considered in his sample (NGC 1399, 
M87, M49) can be accounted for if the mass in GCs is 
normalized to the total baryonic mass of the host (stel- 
lar plus the hot X-ray emitting gas), implying a constant 



baryonic GC formation efficiency e = 0.26%. The parame- 
ter e represents the efficiency of converting baryonic mat- 
ter into GCs, e = M.GCs/Mb- However, as suggested by 
Blakeslee et al. (1997) and Blakeslee (1999), perhaps more 
fundamental is the ratio between the mass in GCs and 
the total galaxy halo mass Mh, (i.e., rjh = Mgcs/ Mh) 
for which they obtain r/ h ~ 1.71 x 10~ 4 . More recently, 
Kravtsov & Gnedin (2005) in a high-resolution cosmo- 
logical simulation for galaxies with M h > IVPMq find 
Mgcs = 3.2 x lO 6 Al0(Xh/lO 11 X0) 13 , which predicts 
en = (2 — 5) x 10~ 5 . Assuming a cosmological baryon frac- 
tion fb — 0.17 (Hinshaw et al. 2009) brings the observed and 
theoretically predicted values of the GC formation efficiency 
to a common value of ~ 10 _J (see Sect. 4.3). Following up 
on the idea of the proportionality between the mass in GCs 
and the host halo mass using a statistical stellar-halo mass 
relation from ACDM simulations, recent studies with much 
larger galaxy samples find a good approximation to the data 
(Peng et al. 2008) and a similar result for the GC formation 
efficiency = 7 x 10~ 5 , Spitler & Forbes 2009). In the 
low-mass galaxy regime, Peng et al. (2008) find evidence for 
an environmental bias: the majority of dwarf galaxies with 
relatively high Sn at a given total luminosity lie within a 
projected radius of ~ 1 Mpc from M 87, the central giant 
elliptical in the Virgo galaxy cluster. The study of Miller & 
Lotz (2007) finds a consistent trend of increasing GC mass 
fraction with decreasing host galaxy mass for Virgo dwarf 
elliptical galaxies. Using the McLaughlin (1999) correction 
for M gas* /M s?ar, Miller & Lotz were able to match the trend 
seen in their observations. 



Investigating the increasing Sn with decreasing galaxy 
mass below M* < 3 x 10 10 M Q , Forbes (2005) used the 
feedback models of Dekel & Woo (2003) which predict 
M/L tx M~' 2 ^ 3 , and found qualitative agreement with the 
observations. However, the normalization of this relation 
remained unconstrained, so that actual GC formation effi- 
ciency values eh could not be determined. This quantitative 
derivation of GC system scaling relations as a function of 
galaxy mass is one of the goals of this work by including 
predictions from the latest Dekel & Birnboim (2006) mod- 
els which include shock-heating regulated star formation for 
M* > 3 x lO 10 ^©. In Section2 we briefly introduce the 
galaxy sample used in this study and access contamination. 
The analysis of the fractions of GCs in late-type dwarf galax- 
ies along with complementing data from the literature is 
presented in Section 3. In Section 4 we discuss the observed 
trends of the various GCS scaling relations, the specific GC 
formation efficiency rj, and implications for galaxy formation 
scenarios. Combined with dynamical mass measurements of 
massive and dwarf galaxies from the literature we normal- 
ize the models (Sect. 4.2.1) and find a good description of 
the observed behavior of the frequencies of GCs as a func- 
tion of galaxy luminosity with a common value of the GC 
formation efficiency parameter. Using the derived analytical 
expressions describing the behavior of the specific GC sys- 
tem scaling parameters as a function of galaxy luminosity we 
discuss (in Sect. 4.4) the predictions of a simplistic satellite 
galaxy accretion model as a function of galaxy luminosity. 
Our final conclusions are summarized in Sect. 5. 
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Figure 1. Contamination analysis. Panels a) and b) show a comparison between the apparent luminosity and color distributions for 
GC candidates in our sample (all, blue and red GC samples with dotted, solid and dashed histograms) and those for sources in the 
GOODS survey (filled histograms) selected similarly to our initial GC candidates selection. Applying identical criteria for the selection 
of contaminants from the GOODS catalog was not possible due to the lack of F814W-band imaging in the latter. Magnitude and 
object morphological selection was based on the GOODS F606W-band measurements and an F606W— F850LP color cut derived from a 
comparison with SSP models, to include the entire range in metallicity at ages > lGyr, matching the F606W— F814W color selection. 
Transformation from F606W— F850LP to F606W— F814W color requires knowledge of the object metallicity. This is well seen in the 
panel c) where we show the F606W-F814W versus F814W-F850LP Galev SSP model tracks (Anders & Fritze-v. Alvensleben 2003) 
for a range of metallicities indicated in the legend. For an approximate comparison, we applied a shift of 0.25 mag to the GOODS 
F606W-F850LP color distribution in panel b). 



2 DATA 

2.1 HST/ACS imaging data 

The analysis of the specific frequencies of mostly late- 
type dwarf galaxies in this study is based on HST/ACS 
archival data from programs SNAP-9771 and GO- 10235 (PI: 
I. Karachentsev), and GO-10210 (PI: B. Tully). These snap- 
shot surveys provide Tip-of-the-Red-Giant-Branch distances 
to all dwarf galaxies in the studied sample (Tully et al. 2006; 
Karachentsev et al. 2006, 2007), and deep observations to 
study the GC systems of the dwarf galaxies. Our sample of 
dwarf galaxies consists of 55 dlrrs, 3 dEs, 5 dSphs, and 5 
Sm galaxies in the field environment, either in loose associ- 
ations of dwarfs or in remote halo regions of nearby groups 
(Sculptor, Maffei 1&2, IC 342, M 81, CVn I cloud). We 
have detected old GCs in 38 dwarf galaxies (cf. Table 1 in 
Georgiev et al. 2009a). Most of the dwarfs have apparent 
diameters smaller than the HST/ACS field of view which 
provides us with a good sampling of their GCSs. The GC 
numbers of two dwarf Sm galaxies from GO-10210 (NGC 784 
and ESO 154-023) and three Sm galaxies from SNAP-9771 
and GO-10235 (ESO 274-01, NGC 247, and NGC 4605) have 
been corrected for spatial incompleteness. Details on the 
data reduction, cluster selection and photometry, complete- 
ness and total galaxy luminosity are described in detail in 
Georgiev et al. (2008, 2009a). 



2.2 Contamination Analysis with HUDF and 
GOODS Data 

An important aspect which can affect the statistics of our 
analysis, especially for the faintest dwarf galaxies with fewer 
detected GC candidates, is the issue of background contam- 
ination. In Georgiev et al. (2008) we used the ACS Hubble 
Ultra Deep Field (HUDF) to select objects with the same 
magnitude and morphological selection criteria as our sam- 



ple GCs. Our analysis showed that the depth and resolu- 
tion of the ACS images allowed us to exclude background 
contamination among the blue GC candidates, even at the 
previously adopted redder cut between blue and red ob- 
jects at (V — I) = 1.1 mag. In our latest studies we have 
adopted a bluer cutoff between blue and red GC candidates 
at (V — J)o = 1.0 mag which is the selection for the cur- 
rent study as well. For a more conservative analysis and 
for further narrowing the possibility of contamination we 
selected for this study GC candidates in the color range 
0.7 < V— I < 1.4 mag, in accordance with the colors of Galac- 
tic and old LMC GCs. 

In order to double-check the contamination level and 
avoid selection biases due to cosmic variance, we have ac- 
cessed the latest version v2.0 of the GOODS catalog (Great 
Observatories Origins Deep Survey, Giavalisco et al. 2004). 
GOODS is a deep imaging survey with HST/ACS in four 
bands (Bf43SW, Vp606W> »F77SW, 2f85ulp) of two locations in 
the north and south sky covering a total of ~320arcmin 2 . 
We matched all sources with photometry in F606W and 
F850LP with a tolerance of 1 pixel for both GOODS-N and 
GOODS-S catalogs. This yields a total number of 39 432 and 
33 955 sources, respectively. 

We analyze the contamination in a similar manner as 
in Georgiev et al. (2008) by employing morphological in- 
formation from the GOODS catalog and selecting sources 
with 19 < F606W < 26 mag, 2 < FWHM F6 o6W < 9 and el- 
lipticity in the range — 0.15, identical to our F606W im- 
ages. By applying the same objects selection strategy as for 
our science images on the entire high-level science GOODS 
v2.0 imaging data 1 we confirmed the morphological selec- 
tion. To constrain the selection of GOODS sources in a sim- 
ilar V—I color range, as the selection of our GC candidates, 



1 Retrieved from http://archive.stsci.edu/prepds/goods/ and 
piped through the same analysis routines of F606W images. 



4 /. Y. Georgiev et al. 



we used the F606W— F850LP color. Observations of other 
galaxies, Galactic and Magellanic Cloud Clusters have John- 
son/Cousins colors in the range 0.7 < V — I < 1.4 mag, which 
we convert using the Galev SSP models (Anders & Fritze- 
v. Alvensleben 2003) to 0.74 < F606W-F850LP < 1.55 for 
ages older than > 1 Gyr for the entire metallicity range. As a 
result, 18 and 24 objects passed the above selection criteria 
for the North and South GOODS fields, respectively. Com- 
bined and normalized to the total GOODS survey area of 
320 arcmin 2 leads to an upper limit of the total contamina- 
tion density of p~0.13 objects per arcmin 2 , i.e. ~ 1.5 objects 
in a single ACS field with an area of 11.33 arcmin 2 . In total 
we have 30 galaxies with blue GC detections (see Sect. 3.1), 
i.e. ~ 340 arcmin 2 which is comparable to the total area of 
the GOODS fields. Thus, virtually no scaling to the number 
of contaminating objects is required for the direct compar- 
ison shown in Figure 1. The luminosity and color distribu- 
tions of all contaminating sources are shown in panels a) and 
b). The magnitude difference between F814W and F850LP 
shown in panel c) is in the range 0.15— 0.45 mag. This shows 
that without knowing the object metallicity it is not possi- 
ble to derive precise information about the contamination as 
a function of color. To partially account for this difference 
and ease the comparison we shifted the color distribution of 
GOODS objects in Figure lb) by 0.25 mag. The luminosity 
distribution of background objects reveals that the contam- 
ination in our entire GCs sample is of the order of ~ 5 — 6 
among the faintest, barely resolved GC candidates. This ad- 
ditional analysis of two different sky regions with deep ACS 
imaging confirms the expected low level of contamination. 
Section 3.1.1 discusses the implications of this result. 



3 ANALYSIS 

3.1 Specific Frequency 

The globular cluster specific frequency is defined as 

S N = A GC x io°- 4 ( M v+is) (1) 

(Harris & van den Bergh 1981) and is a measure of the 
ratio of the formation efficiencies of star clusters to field 
stars. It describes how likely it is for a given galaxy to have 
formed star clusters that survived the various cluster dis- 
ruption mechanisms until today. Two factors affect the Sn- 
values in late-type dlrr galaxies: the recent star formation 
activity of the host galaxy affects My while the dynamical 
evolution history of its GC system affects Agc- 

In particular, for our sample galaxies the influence of 
changes in galaxy luminosity is more significant for late-type 
dwarf galaxies for which the present-day starburst activity 
can boost their luminosity, causing a decrease in the Sn 
values. Formation and subsequent dynamical evolution of 
GCs in low-mass galaxies (where fewer clusters are formed) 
leads to stochastic fluctuations in Agc which can lead to 
an increased dispersion of Sn among such dwarfs, which is 
indeed observed. In addition, the stochastic nature of the 
field-star formation history for dlrrs (e.g. Weisz et al. 2008) 
can introduce additional scatter. 

In spite of this, the field star and star-cluster formation 
efficiencies generally appear to be closely correlated as it is 
observed that the ratio of cluster to star formation rates 
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Figure 2. Globular cluster specific frequency, Sn, versus abso- 
lute galaxy magnitude, My, for all dwarf galaxies in our sample. 
Symbols are the bGC+rGC estimates of the Sn values for vari- 
ous galaxy types, while plus signs indicate the values based on the 
bGC population only.Grey triangles at the bottom of the plot rep- 
resent galaxies in which no GCCs were detected. The black arrow 
indicates the change of the Sn value due to passive evolutionary 
fading of the galaxy light by AMy = 1.5 mag (from 3 to 14 Gyr 
at LMC metallicity). 

(SFR) is nearly constant (e.g. Lada & Lada 2003; Lamers & 
Gieles 2008; Gieles & Bastian 2008) and independent of the 
present-day SFR (Bastian 2008). In addition, the constant 
ratio between star-cluster mass and the mass in field stars at 
the end of violent relaxation 2 (i.e. the bound cluster-to-field 
stellar mass ratio) may be evidence for a nearly universal 
distribution function of the local star formation efficiency 
(Parmentier & Fritze 2009). 

3.1.1 Specific frequencies of dwarf galaxies in our sample 

We compute two Sn values for our sample of GC candidates 
selected from Georgiev et al. (2008, 2009a) in the color range 
0.7 < V—I < 1.4 mag. First using the number of blue GCs and 
second the total number of GCs: blue GCs + red GCs (red 
GCs have (V—I)o > 1.0 mag), together with the total V-band 
galaxy magnitudes derived in Georgiev et al. (2008, 2009a). 
The calculated Sn values for all dwarf galaxies in our sample 
are shown in Figure 2 and listed in Table 1. The contamina- 
tion analysis in Section 2.2 suggests that red GC candidates 
with 1.0 < (V—I) < 1.4 can not be fully accounted for by con- 
tamination in our sample galaxies. This is puzzling because 
the existence of such red, metal-rich GCs in low-mass dwarf 
galaxies is not expected. The fraction of red GCs in early- 
type galaxies decreases with decreasing galaxy luminosity 
and becomes negligible for the faintest dEs (cf Fig. 4 and 10 
in Peng et al. 2006). Theory suggests that red, metal-rich 
clusters are believed to form either in major galaxy mergers 
(e.g. Ashman & Zepf 1992) or in a multi-phase collapse of 

2 The end of violent relaxation, when the dynamical response of 
a cluster to the expulsion of its leftover star forming gas is com- 
pleted, typically takes place at an age between 10 and 50 Myr and 
depends on the external tidal field and the proto-cluster crossing- 
time (see Fig. 4 in Parmentier 2009). 
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a massive galaxy (Forbes et al. 1997; Pipino et al. 2007). 
Both channels are unlikely to occur in low-mass systems. 
Observationally, extinction within the galaxy can lead to a 
redder colour. Taking the Magellanic Clouds as a proxy for 
the dwarfs in our sample, the expected internal reddening is 
in the range < E(V-I) < 0.13 mag 3 The radial distribu- 
tion of the bGCs and rGCs in our sample galaxies are indis- 
tinguishable from each other. They are distributed mostly 
within 1-2 pc projected distance from the galaxy center, i.e. 
~ lxR c ff (Georgiev et al. 2009b). Therefore, some inter- 
nal reddening can be expected among the GCs in our sam- 
ple. Red colors could arise if clusters have ages in the range 
1 — 3 Gyr where AGB stars are in their thermal-pulsation 
phase (TP-phase) which affects the integrated cluster color. 
Comparison with population synthesis models (e.g. Maras- 
ton 2005; Leitherer et al. 1999) shows that for that age range 
the integrated cluster colour can reach V — I < 1.4 mag for 
even as low metallicities as 0.2 to 0.4 Z@. However, whether 
reddening or younger ages are affecting the colors of GCs in 
our sample, or whether they are background contaminants 
will require spectroscopic or near-IR follow-up observations 
to establish their accurate metallicities and membership. 

In order to minimize the effects on the Sn from the un- 
certain nature of the rGCs in our subsequent analysis, we 
consider only dwarf galaxies which contain blue GC can- 
didates. There are eight galaxies (20%) with only red GCs 
which are treated as galaxies with null GC detections. The 
lack of bGCs does not necessarily mean that some of the 
rGCs cannot be real GCs due to the effects discussed above. 
Therefore, given these uncertainties, our Sjvs represent a 
conservative lower value for faint late-type galaxies. In Sec- 
tion 4.2.2 we discuss the implications and limitations of the 
presence of rGCs for the model we are testing. 

Overall, the observed trend of increasing Sn with de- 
creasing host galaxy luminosity holds for the entire sample, 
and in particular for dlrrs. Because of their small sample size 
the trends for dEs/dSphs and Sm galaxies are less conclusive 
from our data alone. Grey triangles show galaxies where no 
GC candidates are detected. Their Sn values were adjusted 
for illustrative purposes only. Notably, Figure 2 shows that 
for low-luminosity galaxies the difference between total Sn 
(where the rGCs sample can include background contami- 
nants) and blue GC Sn (the most trustworthy value) differ 
a little quantitatively, but not significantly in the qualita- 
tive sense due to the logarithmic scaling. Thus, the observed 
general trend of increasing Sn value with decreasing galaxy 
luminosity is present even if only Sn values calculated from 
the blue GCs (bGCs) are considered. 

It is important to note that all Sn values for the late- 
type dwarfs represent their present-day values. Passive evo- 



3 The total E(B — V) value for old LMC GCs is in the range 
0.083 - 0.213 mag for NGC 1898 and NGC1841 (McLaughlin & 
van der Marel 2005). Corrected for the median Galactic fore- 
ground reddening toward LMC of E(B— V) = 0.075 mag (Schlegel 
et al. 1998), leads to E(B-V) = 0.008 and 0.138mag, i.e. 
E(V-I) = 0.009 and 0.164 mag internal reddening for the LMC. 
The same analysis for the only old SMC GC, NGC 121, gives 
E(B—V) = 0.109 mag or E(V-I) = 0.130 mag. The LMC GCs 
above are at R = 0.6° and 14.9° projected distance from the cen- 
ter of the LMC (van den Bergh & Mackey 2004). m-M = 18.5 mag 
(50.118 kpc) to the LMC yields 0.5 and 13.3 kpc. 
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Figure 3. GC specific frequency as a function of absolute galaxy 
magnitude for dwarf galaxies at low density environment from 
Sharina et al. (2005); Georgiev et al. (2008, 2009a) and dlrrs in 
Virgo and Fornax cluster Seth et al. (2004); Georgiev et al. (2006). 
Dwarf irregular galaxies are shown with solid circles, plus signs 
and asterisks; dSphs with open and dashed-open circles; Sms with 
solid squares. With solid pentagons we show Sn values of the best 
studied GC systems in nearby dlrrs — the LMC and SMC. Grey 
triangles at the bottom of the plot represent galaxies in which no 
GCCs were detected. The arrow indicates passive evolutionary 
fading as in Figure 2. 

lutionary fading of the integrated galaxy light will increase 
the Sn values by a factor of 2 to 16 for our dlrr and 
Sm galaxy sample (assuming Nqc = const.). This is shown 
by an arrow in Figure 2, which indicates the change in 
Sn if a galaxy passively fades due to stellar evolution by 
AMy = 1.5 mag. This corresponds to an evolution from 3 
to 14 Gyr at LMC metallicity as inferred from SSP models 
(e.g. Bruzual & Chariot 2003). 



3.1.2 Comparison of Sn values of late-type dwarf galaxies 

In Figure 3 we compare the Sn values of our sample dwarfs 
with those of late-type dwarfs from Sharina et al. (2005), 
all in field environment, and Virgo and Fornax cluster dlrrs 
from Seth et al. (2004) and Georgiev et al. (2006), as well as 
the Magellanic Clouds. It is clearly seen that, combined with 
the data from the literature, we observe a general trend of 
increasing Sn with decreasing galaxy luminosity for dlrrs. 
Such a trend has been previously observed for early-type dEs 
(e.g. Miller et al. 1998; Strader et al. 2006). The dlrrs in our 
sample are systematically brighter than dSphs by ~ 3 mag, 
indicating younger ages. As mentioned above, a simple pic- 
ture of passive evolutionary fading of the dlrrs can shift their 
Sn values toward the dSphs region, parallel to the average 
trend. This supports the notion that if the number of GCs 
remains unchanged the high— z analogs of present-day dlrrs 
might be the progenitors of the more evolved dSphs/dEs 
(e.g. Miller et al. 1998; Seth et al. 2004). 

To investigate the formation efficiency of GCs as a func- 
tion of galaxy luminosity, we extend the galaxy mass range 
and morphology by including data from the literature. In 
Figure 4 we compare the specific frequencies of the dwarf 
galaxies discussed above with those of early-type galaxies in 
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Figure 4. Top: GC specific frequency versus absolute galaxy 
magnitude for a range of galaxy masses, morphologies, and envi- 
ronments. The various symbol types are explained in the figure 
legend. Massive and dwarf elliptical galaxies as well as spirals 
have been collected from the literature as explained in Sect. 3.1. 
With larger symbol sizes (dots, circles and squares) we indicate 
the dwarf galaxies in our sample. Grey triangles at the bottom of 
the plot represent galaxies in which no GCCs were detected. The 
large symbols connected with a line are the co-added running av- 
erage of the Sn values per magnitude bin for the different galaxy 
morphological types. The errorbars are the standard deviation at 
that magnitude bin. The dash-dotted line indicates the Sn value 
if a galaxy hosts one GC (cf. Eqn. 5). Bottom: Fraction of galax- 
ies with GC detections for the magnitude bins used to derive the 
mean Sn value in the upper panel. 



the Virgo galaxy cluster from Peng et al. (2008) and other 
massive elliptical and spiral galaxies from the compilation by 
Spitler et al. (2008) who selected galaxies with robustly de- 
rived GCS properties (deep photometry and completeness 
limits, adequate spatial coverage etc.) from the literature 
(Ashman & Zepf 1998; Rhode & Zepf 2001, 2003, 2004; 
Forbes et al. 2001; Forte et al. 2001; Dirsch et al. 2003, 2005; 
Gomez & Richtler 2004; Bassino et al. 2006b, a, 2008; Harris 
et al. 2006; Tamura et al. 2006a,b; Rhode et al. 2007). Our 
main source for the early-type dwarf galaxy data is Miller & 
Lotz (2007). For reference, we have labeled in Fig. 4 the cor- 
responding Sn values of the Magellanic Clouds, M31, and 
the Milky Way. Keeping in mind the effects of passive evo- 
lutionary fading and star cluster disruption, the Sn vs. My 
trends for late- and early-type dwarfs in cluster and in loose 
environments appear very similar irrespective of galaxy mor- 
phological type and environment. With a dash-dotted line 
we indicate the lower limit of the Sn value defined by Equa- 
tion 1 if the galaxy hosts one GC (Ngc = 1)- 

To investigate the Sn value as a function of galaxy mor- 
phology and luminosity we have estimated the mean Sn 
value in bins of one magnitude. With large open symbols 
connected with lines in Figure 4 we show these co-added 



running averages of the Sn including galaxies without GC 
candidates. They have been computed by summing the num- 
ber of GCs and galaxy luminosity for each magnitude bin. 
The error bars represent the standard deviation at that mag- 
nitude bin. It can be seen that late-type galaxies (spirals and 
irregulars) have on average (~ 2 times) lower mean Sn val- 
ues at a given galaxy luminosity than early-type galaxies 
(ellipticals and spheroidals). This constant offset can be ac- 
counted for by evolutionary fading effects. The spiral galax- 
ies in our sample show a fairly constant Sn value as a func- 
tion of galaxy luminosity (Sn ~ 1). However, a larger sample 
increasing the number of spiral galaxies with luminosities 
Mv < —19 mag need to be studied to confirm the nearly 
constant Sn— value for spiral galaxies. 

The bottom panel of Figure 4 shows an increasing frac- 
tion of galaxies with no detected old GCs towards decreas- 
ing galaxy luminosities. This decrease can be attributed to 
several factors affecting the galaxy luminosity normalized 
number of GCs: i) The increased impact of stochasticity for 
star-cluster formation and/or ii) the decrease of GC num- 
bers due to dynamical destruction processes (i.e. evapora- 
tion, tidal disruption, etc.). 



3.2 Specific Mass, Luminosity, and T-value 

In addition to the traditionally used specific frequency Sn, 
we discuss in the following the trends in the specific mass and 
specific luminosity of the GC system, as well as a redefined 
T-value (see Zepf & Ashman 1993) as a function of host 
galaxy luminosity and mass. 



3.2.1 Specific Mass of GC system 

We define the "specific mass" Sm as the total mass of the 
globular cluster system, A(gcs, relative to the total baryon 
mass of the host galaxy defined as the sum of the stellar 
mass, Al*, and the HI gas mass, Aim: 



Sm = 100 x Mgcs/{M* + Mm). 



(2) 



This is similar to 100 times the McLaughlin (1999) e param- 
eter for the GC formation efficiency for giant galaxies with 
respect to the total baryon mass, which includes the mass 
of the hot, X-ray emitting halo. It is not a trivial task to 
include this mass component for the galaxies in our sample. 
Obtaining hot gas mass from X-ray profiles (Mx) relies on 
assumptions of the temperature profile as well as hydrostatic 
equilibrium, which is not always the case for galaxies in com- 
plex environments. For instance, the Fukazawa et al. (2006) 
X-ray observations of luminous early-type ellipticals in a 
range of environments show much weaker X-ray emission 
from isolated than for clustered galaxies. In addition, they 
find that the requirement for hydrostatic equilibrium might 
not be fulfilled for many of their ellipticals due to jets from 
AGNs and interaction with neighboring galaxies. In general, 
all cD galaxies in their sample (like NGC 1399) are classified 
as galaxies embedded in the hot cluster/group ISM. There- 
fore, the X-ray halo of such central cluster galaxies may have 
a non-negligible contribution from the galaxy cluster itself. 
To avoid problems with hydrostatic equilibrium Humphrey 
et al. (2006) and O'Sullivan et al. (2007) studied isolated el- 
liptical galaxies for which they obtain Mx few times 10 8 Mq, 



a similar Mx mass was estimated for M49 by McLaughlin 
(1999). Therefore, for the majority of the galaxies in our 
sample the expected Mx is of the order of few percent of 
the stellar mass, thus a negligible contribution to the total 
galaxy mass. 

In the following we have considered the sum of the HI 
and stellar masses, M*+Mhi, for the dwarf galaxies in our 
sample. Including the HI gas mass becomes more important 
when comparing massive elliptical galaxies with dwarf galax- 
ies which can contain significant neutral gas mass fractions. 
The disadvantage of the Sm parameter (Eq. 2) is the neces- 
sary conversion from luminosity to mass, using M/L ratio 
for the galaxy as well as the GCs (see below). On the other 
hand, Sm has the advantage of being a distance-independent 
quantity. 

To derive .Mgcs for dlrrs we convert GC absolute mag- 
nitudes to masses using 7v = (M/Lv)gc = 1.88, which is the 
mean V-band mass-to-light ratio for old GCs in the LMC 
that we estimated from McLaughlin & van der Marel (2005). 
Since we are studying old globular clusters in similar galax- 
ies, our choice of 7v is justified and changes only within 
10% for the age range of our GC sample selection (t > 4 
Gyr). However, a recent study by Kruijssen & Portegies 
Zwart (2009) showed that the GC M/L ratio is expected to 
depend on cluster luminosity (mass) in that low-mass clus- 
ters preferentially lose low-mass, i.e. high M/L, stars due to 
two-body relaxation. They show that this effect influences 
the translation from GC luminosity to mass function and 
the respective location of the turn-over mass (m-ro) toward 
lower mass by ~ 0.2 dex. This effect needs to be tested for a 
range of galaxy masses (i.e., tidal field strengths) in order to 
establish the GC M/L oc L effect to the GCS m/po- Keep- 
ing this in mind and assuming a constant GC M/L for the 
GCSs of all galaxies, we will show in Sect. 4 that this effect 
is negligible relative to the observed spread/scatter in the 
relations that involve the GC M/L ratio. 

Mm was calculated from the HI magnitudes listed in 
the LEDA database (Paturel et al. 2003) or from HI fluxes 
in NED or Begum et al. (2008) following Roberts & Haynes 
(1994), i.e. Mm = 2.36 x 10 5 x D 2 x F m , where D is the 
distance in Mpc and Fm is the HI flux in Jy kms" 1 . To 
convert from galaxy luminosity to galaxy stellar mass M*, 
we used the relations in Bell et al. (2003) between the galaxy 
B — V colour and the M+/Lv ratio (or B — R if B — V was 
not available): 

log(M*/L v ) = -0.628+ 1.305 x (B-V) 
\og{M*/L v ) = -0.633 + 0.816 x (B-R). 

In addition to our dlrr galaxy data, we calculated Mi,/ Ly 
ratios for elliptical and spiral galaxies from Peng et al. 
(2008); Spitler et al. (2008); Miller & Lotz (2007) and for 
dSphs and dlrrs from Sharina et al. (2005) and Seth et al. 
(2004) as follows. We use the stellar masses directly from the 
Peng et al. study which were computed from galaxy colours 
(g — z and J—K s ). Spitler et al. (2008) used A'— band pho- 
tometry and theoretical M*/Lk from population synthesis 
models to derive galaxy mass, which we also adopt in the 
subsequent analysis. To derive baryonic masses for the dE 
sample in Miller & Lotz (2007) we use the galaxy magni- 
tudes provided in their table and M+/Ly = 3.0. We used 
this mass-to-light ratio because the masses listed in their 
table were calculated assuming M*/Lv = 5 which seems 




-22 -20 -18 -16 -14 -12 -10 -8 
M v 

Figure 5. Stellar plus HI mass as a function of galaxy luminosity 
for dwarf galaxies in our sample (large solid dots) compared with 
galaxies of different morphological type from the literature as 
indicated in the legend. Smaller solid dots show dlrrs from Seth 
et al. (2004) and Sharina et al. (2005). With solid pentagons are 
shown (frm left to right) M31, Milky Way, LMC and SMC. The 
dashed line is a least-squares fit to the data; the fit relation and 
the values of the coefficients are given in the lower left corner. 
The solid line shows the relation between the galaxy luminosity 
and its mass with M t /Ly = 3.0. 



to be too high and offset from the relation defined by our 
dataset and those of Seth et al. (2004); Sharina et al. (2005); 
Peng et al. (2008); Spitler et al. (2008) (see Fig. 5). For the 
dwarf galaxies in the Sharina et al. (2005) and Seth et al. 
(2004) samples we calculate their masses from the Bell et al. 
relations as we did for our sample. Note that uncertainties in 
the galaxy M/L determination due to unknown SFHs and 
the choice of an IMF apply equally to the range of galaxy 
masses in our sample as to that of the SDSS and 2MASS 
galaxy samples from which the Bell et al. relations were de- 
rived. 

In Figure 5 we illustrate the relation between galaxy 
luminosity and its stellar plus HI mass, M*+m, for galaxies 
from the literature and the dwarf galaxies in our sample. 
We find an increasing difference with decreasing luminosity 
between the late-type, gas-rich, low-mass dlrrs, and early- 
type dwarf ellipticals, which are essentially devoid of HI. 
In addition, the low-mass dlrrs show an increasing disper- 
sion towards fainter systems. This offset and dispersion can 
both be attributed to possibly underestimated Mb/Ly due 
to younger stellar populations and/or unaccounted-for bary- 
onic components in the form of, e.g. diffuse ionized gas, hot 
X-ray emitting gas, radio continuum emission or molecular 
gas. All these mass components are known to be associated 
with star formation activity and their mass fractions are 
observed to scale with galaxy HI mass (e.g. Dahlem et al. 
2006). 

3.2.2 Specific Luminosity of GC system 

Another useful indicator of GC formation efficiency is the 
specific GC luminosity introduced by Harris (1991) 

S L ,v = 100 x L G cs/Lv (3) 
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which is the luminosity ratio between the GCS, Lgcs, and 
its host galaxy, Lv- Its main advantage as scaling parameter 
is that it is independent of distance and only weakly sensi- 
tive to completeness corrections at the faint end of the GC 
luminosity function (GCLF). A slight disadvantage of Sl is 
its sensitivity to the stellar population dominating the inte- 
grated light of the galaxy. For the case of late-type galaxies 
this typically is the most prominent (the latest) star forma- 
tion episode, so that the galaxy luminosity does not neces- 
sarily indicates the stellar population at the time when most 
of the GCs were formed. On the other hand, Sl is indepen- 
dent of M / L and hence least affected by stellar population 
systematics. Thus, Sl may be the most robust GCS scaling 
parameter for comparing GCS richness and GC formation 
efficiencies between galaxies over large range of masses. 

3.2.3 Modified T -value 

Zepf & Ashman (1993) introduced the T value as the spe- 
cific number of globular clusters per unit 10 9 Mq of the host 
galaxy mass, i.e. T = 1O 9 M x N G c/M*. For our purposes 
this definition is insufficient since we want to compare gas- 
depleted early-type with late-type galaxies with significant 
gas mass fractions. Thus, we modify the T- value from Zepf & 
Ashman by including the gas mass of the host galaxy. Hence 
we define the specific number of GCs per unit baryonic mass 
f (Mb) as follows: 

f = 10 9 M Q x Nac/Mt (4) 

where Mb = M*+Mm- 

In Figure 6 we directly compare the GC specific fre- 
quencies (Sn, Eq. 1), the mass-normalized specific mass in 
clusters (Sm, Eq. 2), specific luminosities (Sl, Eq. 3), and 
the specific number (T, Eq. 4) as a function of host galaxy 
luminosity for the entire sample. 



4 DISCUSSION 

In this section we address the question whether there is a 
physically motivated model which is able to describe the ob- 
served trends of the various scaling relations as a function 
of total galaxy mass. This can be tested by assuming that 
GCs form in proportion to the total galaxy mass. We start 
by comparing the GC scaling relations introduced in the pre- 
vious Section 3 with theoretical model predictions that pro- 
vide relations between the mass-to-light ratio and the total 
galaxy mass (Dekel & Silk 1986; Dekel & Birnboim 2006). 
The star-formation processes in these models are driven by 
(i) stellar feedback and (ii) shock heating below and above 
a critical stellar mass of 3 x 10 lo .M©, respectively. Thus, 
we derive analytical expressions that can be used to con- 
front the observed trends against the assumptions of con- 
stant and variable GC formation efficiencies as a function of 
galaxy luminosity. Observed differences and systematics can 
argue for or against physical mechanisms included (or the 
lack thereof) in the models we discuss. 

Is the assumption of an universal mode of globular clus- 
ter formation among galaxies of different morphological type 
and environment justified? An universal GC formation sce- 
nario would predict that all (evolved) galaxies form their old- 
est star cluster populations with the same initial efficiency 



at given galaxy mass. Subsequent formation of metal-rich 
GCs in massive galaxies would likely cause a spread in the 
observed GC formation efficiency value for a given range in 
galaxy mass since galaxies undergo different merging his- 
tories, varying from predominantly dissipational to mostly 
dissipationless. This can be tested by comparison of the ex- 
pected value of the GC efficiency from the adopted analyti- 
cal models and the observations of galaxies of the two major 
morphological types, early- versus late-types, and as a func- 
tion of galaxy mass. The relative difference between the two 
galaxy types should be preserved at a given galaxy mass, 
assuming that GC destruction processes are similar to first 
order in both galaxy types. 

In order to set quantitative constraints and limitations 
of the models we use, and address all questions discussed 
above, it is very important to normalize the model predic- 
tions, i.e. to put the models on an absolute scale. All studies, 
so far, which addressed the shape of the Sn distribution used 
one or another version of the Dekel & Birnboim models (e.g. 
Forbes 2005) but did not derive strict constraints due to the 
lack of such a normalization. Therefore, our effort in deriv- 
ing a normalization will help to narrow down the possible 
physical processes that govern the spread in GC formation 
efficiencies among galaxies with similar mass. This is de- 
tailed in Section 4.2.1. 



4.1 Boundaries of Globular Cluster Systems 
Scaling Relations 

The definition of Sn (Eq. 1) determines the most trivial 
lower limit of star cluster population when Ngc = 1, i.e. 
the GC formation efficiency that resulted in only one old 
GC at present (and no additional field star component), so 
that 

logSV = 6 + OAMv = 7.928 - log L v , (5) 

assuming Mv,@ = 4.82 mag (Cox 2000). Analogously, we 
can derive similar relations for Sl,Sm and T by assuming 
that the total luminosity and mass of the globular cluster 
system are Lgcs = A^gc x ho and Macs = Nqc x m T o, 
where ho and mxo are the cluster luminosity and cluster 
mass at the turn-over of the GCLF, respectively, hence 

log Sl = 2 + log /to — log L v 

= 2 + log ho + 0.4 (M v - M v ,@) (6) 

log 5m = 2 + log(m T o/Tv) - logLv 

= 2 + log(m T o/Tv) + O.4(Mv-M v , ), (7) 

log T — 9 — log T v ~ log L v 

= 9-logTv + O.4(Mv-M v , ) (8) 

where Tv = M*/Lv is the (overall) T/-band stellar mass- 
to-light ratio of the host galaxy. The relations (5) to (8) are 
shown in Figure 6 with dash-dotted lines. For the derivation 
of Sm and f we have adopted T v = 3.0, ho = 8.47x 10 4 L Q , 
and m.TO = 1-69 x 1O 5 M0, which corresponds to Mv = 
—7.5 mag for 7y = 2.0. Note that a change from Tv = 3.0 to 
1.0 and from Mv = —7.5 to —7.4 mag causes a parallel shift 
comparable to the symbol size relative to the relations. 
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Since the quantities Sn, Sl, Sm, and t involve prop- 
erties of GC systems that are (often) measured down to the 
turnover luminosity of the GC luminosity function, it is use- 
ful to assess the uncertainties associated with Eqs. 5-8. The 
results of Jordan et al. (2007) show that the dispersion (a) 
and turn-over luminosity (/to) of the GCLF changes as a 
function of host galaxy luminosity. The same study shows 
also that the mass at the GCLF turnover is roughly constant 
(see also Vesperini 1998, 2000) and decreases by ~ 30% from 
rriTO = (2.2±0.4) x 1O 5 M in the brightest galaxies down to 
m-ro — (1.6—1.7) x 1O 5 M in dwarf elliptical galaxies, with 
a substantial scatter in the latter. For comparison, in our 
analysis of the GCLF in mainly late-type dwarf galaxies we 
derive a GCLF turnover mass mTO = 1.6x10 Mq (Georgiev 
et al. 2008) assuming yv = 1.88, which is the mean value 
of the mass-to-light ratio of old GCs in the LMC and SMC 
dlrrs (McLaughlin & van der Marel 2005). Due to the fact 
that the HST/ACS field of view entirely covers the dwarf 
galaxies in our sample we observe the total GC population 
and derive an accurate total Sn 4 - Doing so we find that the 
predictions of the Equations (5) to (8) are accurate to within 
~ 30% for our sample of late-type dwarf galaxies. 

The above relations can be understood as a minimum 
GC formation efficiency threshold. In other words, the in- 
verse scaling relations with galaxy Ly suggest that dwarf 
galaxies at fainter total luminosities are either more efficient 
in forming star clusters relative to the formation of field stars 
or less efficient in the formation of field stars than more lu- 
minous systems (see Fig. 6). We have discussed above (see 
Sect. 3.1) that passive luminosity fading of the host is al- 
most parallel to this lower limit and leads to an increase in 
the Sn value, provided the GC population size remains con- 
stant. However, it should be noted that the slope observed 
in the data partially reflects the slope stemming from the 
definition of these equations. 

4.2 Globular Cluster System Scaling Relations as 
a Function of Galaxy Mass 

All panels in Figure 6 strengthen the previously observed 
increasing GC formation efficiency with decreasing galaxy 
luminosity for galaxies fainter than My « —20 mag. This 
trend is preserved even when the GC number (or mass) is 
compared to the total host galaxy mass. Theoretical models 
of star-formation processes regulated by supernova feedback 
in low-mass halos (e.g. Dekel & Silk 1986; Dekel & Woo 
2003), predict an increasing total Mh/L with decreasing 
halo mass. The assumption that the number of GCs could 
be a function of host galaxy mass-to-light ratio and environ- 
mental density have been investigated in previous studies of 
early-type galaxies (e.g. Forbes 2005; Bekki et al. 2006; Peng 
et al. 2008, and references therein). 

In the following we derive analytical expressions which 
relate the observed specific formation efficiency of old GCs 
in galaxies as a function of total halo mass Mh- Under the 
assumption that GCs form in proportion to the mass of the 

4 Note that the original definition of Harris & van den Bcrgh 
(1981) is to double the number of GCs brighter than the GCLF 
turnover luminosity to obtain A^qc in Equation 1. This is identical 
to the total Nqq only if the GCLF is symmetric 



host galaxy, i.e. the total GC system mass A4gcs is directly 
related to M.h as: 

-Mgcs = nMh, (9) 

where n represents an empirical formation efficiency param- 
eter, which is the specific globular cluster formation effi- 
ciency. This parameter describes the total observed mass 
fraction in globular clusters, which is a product of the ini- 
tial cluster formation efficiency (initial cluster mass func- 
tion) and various mechanisms of cluster dissolution (e.g. 
two-body relaxation, GMC encounters, tidal shocking and 
dynamical friction). Note that our sample dwarf galaxies 
have old GCs with ages expected from their V — I colour 
of t > 4 Gyr, hence, we are probing the early, post-violent 
relaxation phases of GC formation and their dissolution his- 
tory with respect the field stellar population. In the follow- 
ing, we assume that the total GC system mass and its cor- 
responding total luminosity are related through 

Lgcs = Mgcs/jv (10) 

where jy is the mean V-band GC mass-to-light ratio. 

The observed GC system scaling relations have been 
tied by previous studies (e.g. Forbes 2005) into the galaxy 
formation context using the models of Dekel & Silk (1986) 
according to which star formation in dark-matter halos be- 
low a stellar mass of M « 3 x 10 10 M Q (Mh ~ 1O 12 M ) 
is governed by the thermal properties of the inflowing gas 
which are primarily regulated by supernova (SN) feedback. 
These models predict that the mass-to-light ratio evolution 
of mini-halos follows M h /L v ocM,; 2/3 (Dekel & Woo 2003), 
so that LyocM h . 

The latest models on star formation regulated by shock- 
heating (Dekel & Birnboim 2006) predict that above this 
characteristic stellar mass (M « 3 x 1O 1O .M ), the gas in- 
flow experiences a virial shock ("hot stream") and, hence, 
is more susceptible to feedback leading to suppressing star 
formation. This physical transition leads to a different mass- 
to-light ratio scaling that follows Mh/Ly ccM]/ 2 or Ly oc 
M\' 2 for M > 3 x 1O 1O 7W . 

4-2.1 Normalization of the Scaling Relations 

In order to obtain a quantitative description of the specific 
GC formation efficiency as a function of galaxy halo mass 
we must normalize these scaling relations. We, therefore, 
assume the following normalization relations 

mL v =Ml /3 for M<3 x 10 10 M Q (11) 

k 2 L v = M\' 2 for M > 3 x 1O 1O X . (12) 

Such a normalization has not been investigated by previ- 
ous studies (e.g. Forbes 2005). To empirically calibrate the 
model predictions in the low-mass, SN-feedback regulated 
star formation regime (X<3x 1O 1O M ), we use the results 
from Strigari et al. (2008) who measure dynamical masses 
within the inner 0.3 kpc of many Local Group dwarf galax- 
ies over a wide range of luminosities (10 3 — 1O 7 L ). Strigari 
et al. find that all dwarfs have a "universal" central mass of 
A4o.3k pc = 1O 7 M . This central dynamical mass can be tied 
through high-resolution cold dark matter (CDM) simula- 
tions to the total halo mass via Mh — ■Mo / 3kpc x 1O _11 M , 
based on the models of Bullock et al. (2001). Below the 
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Figure 6. Globular cluster system scaling parameters as a function of galaxy luminosity. From top left to bottom right are shown the GC 
specific specific luminosities (Sl), specific mass (Sm), specific frequencies (Sjv), and specific number (T) for all galaxies in the combined 
sample. Large solid dots indicate dlrrs from our study, while small solid dots are dlrrs from Seth et al. (2004); Sharina et al. (2005) and 
Georgiev et al. (2006). Solid pentagons show the corresponding values for the Milky Way, M31, LMC and SMC. Grey triangles show 
galaxies in which no GC candidates are detected. The dash-dotted line represents the corresponding value if galaxy hosts one GC at 
present, see Sect. 4.1. The solid curves are predictions by the models which assumes that GCs form proportional to the total galaxy halo 
mass and that stellar, SNe-driven feedback and virial shock- heating regulate galaxy wide star formation (Dekel & Birnboim 2006), below 
and above M ~ 3 X 10 10 .M©, respectively (see Sect. 4.2). The solid curve (top panels) is the best-fit t)l = 5.5 and t)m = 6.7 X 10 -5 
value to the Sl and Sm distributions, respectively. The solid lines in the bottom panels are curves with adopted t)l = 5.5 X 10 . 
Large triangles connected with solid lines in the left panels show the co-added running average Sjv and Sl values as a function of 
galaxy luminosity for early- and late-type systems. Dashed curves illustrate the range in r) values among galaxies, between r] = 10~ 5 to 
3.5 X 10" 4 . 



critical mass of M = 3 x 1O 1O .M0 this leads to k\ = 
(Xo^kpoXHT 11 ) 573 /Lv- For L v , we choose L V = W 6 L Q as 
the lowest luminosity of a stellar system that can possibly 
form one GC. Furthermore, the dynamical mass measure- 
ments of Strigari et al. are less dominated by observational 
errors at this luminosity. For the high-mass, "virial-shock" 
regime we base the determination of K2 on empirical results 
for the virial mass to B-band light ratio at 3 X 10 10 .M©: 
Mvii/Lb = 75 from Eke et al. (2006) based on galaxy group 
dynamics. From SSP models at solar metallicity (Bruzual 
& Chariot 2003) we derive (M / L B ) / {M / L v ) = 1.65 at 
14Gyr which leads to M v i v /Lv = 41. Hence, we obtain 
L v = 6.67 X 10 8 X©, i.e., « 2 = 10~ 4 42 for the high-mass 
regime. Thus, we have 
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(14) 



Using Equations (9) through (11) we can express the scal- 
ing relation defined by Equations (1), (2), (3), and (4) as a 



function of galaxy absolute magnitude for the "SN feedback" 
regime at stellar mass Mh <3 x 10 10 .M©: 

logSV ~ 6 + 0.24 x M v ,q +log +0.16M v (15) 
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m TO 

( rati* 

V 7v- 



log Sm ^ 2 - 0.16 x M v , e + log 



riK 06 
r l li l 
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+ 0.16M V (17) 



logT ~ 9 - 0.16 x M 



> fl";;| VK °'l ) +0.16My- (18) 



where Tv is the baryonic l/-band mass-to-light ratio of the 
host galaxy, rriTO is the GCLF turn-over GC mass, and jv is 
the corresponding GC mass-to-light ratio. The scaling rela- 
tions as a function of galaxy luminosity Ly, total halo mass 
A4h, and baryonic mass Mb are given in Appendix A. 

For the "virial shock" regime at masses M > 3 x 
lO 10 ^© {Mb, > 10 12 A4©) the above scaling relations have 
the following functional forms: 
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(19) 
(20) 
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4-2.2 Continuous description of the GCS scaling relations 

In order to simultaneously describe the observed trends 
in the various GCS scaling relations in the stellar/SNe- 
dominated feedback regime and the regime where star- 
formation is primarily regulated by virial shock-heating of 
the inflowing gas, we combine the above sets of relations 
above and below the characteristic galaxy stellar mass at 
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Solid and dotted curves in Figure 6 illustrate the rela- 
tions from Eqs. (23) to (26), for various values of r\. It can be 
seen that these scaling relations approximate the observed 
distributions fairly well and indicate a statistically signifi- 
cant spread in the mass and luminosity normalized GC for- 
mation efficiencies t)l and t\m (top panels in Fig. 6). The 
co-added running average Sl and Sn values as a function 
of galaxy luminosity for early- and late-type systems (con- 
nected with solid line large triangles in the left panels of 
Fig. 6) also show that they compare well with the analytical 
model expectation. 

The remarkably good description of the observations by 
the analytical model advocates that at low galaxy masses, 
where galaxies have shallower potential wells (i.e., lower 
binding energy), gas is effectively heated up by the ioniz- 
ing radiation of SNe suppressing star-formation. In the view 
of the Dekel & Birnboim models star-formation in massive 
and dense molecular clouds is effectively shielded from heat- 
ing radiation which does not significantly influence cluster 
formation but suppresses more effectively the formation of 
field stars which leads to an increase in the GCS scaling pa- 
rameters. Thus, at low-galaxy masses the initial burst of star 
formation leads to the formation of the old metal-poor GCs. 
The subsequent cooling of the galaxy ISM and its chemi- 
cal enrichment can lead to the formation of metal-rich and 
younger GCs triggered by instability of internal (tidal den- 
sity waves and/or collisions of stellar winds) or external (e.g. 
galaxy-galaxy interaction) origin. Examples are the metal- 
rich ([Fe/H] < —0.5) and young (>fewGyr) clusters in the 
Magellanic Clouds and Local Group (LG) dwarf galaxies. 



This is consistent with the recent age-metallicity relation 
(AMR) derived for those systems (Forbes & Bridges 2010; 
Sharina et al. 2010). These studies show that the formation 
of star clusters in low-mass systems captures a snapshot of 
the enrichment history of the host galaxy which is a func- 
tion of host mass expressed by the metallicity-luminosity 
relation of dSphs/dlrr galaxies (e.g. Mateo 1998; Crnojevic 
et al. 2010; Lianou et al. 2010). If some of the red GCs in 
our sample are indeed young (see discussion in Sect. 3.1.1) 
this might support this scenario. 

With increasing galaxy mass, the galaxy potential deep- 
ens, the binding energy is larger, which leads to a more ef- 
fective field star-formation efficiency at the same initial GC 
formation efficiency, thus causing the decrease in scaling pa- 
rameters, e.g. Sl- Beyond the virial shock-heating galaxy 
mass, dense molecular clouds are better shielded than the 
lower-mass, less compact gas clouds in which most of the 
field star-formation will occur. On the one hand, this will 
lead to a nearly universal initial GC formation efficiency, on 
the other hand, to a more efficiently suppressed field star- 
formation efficiency. This will cause the observed upturn in 
{23) Sl at high galaxy masses. That is, at the transitional stellar 
mass of M ~ 3 x 1O 1O 7M0 the thermal properties of the 
inflowing gas are primarily governed by virial shocks which 
heat the ISM and inhibit the efficient formation of field stars 
at the same initial GC formation efficiency. The variance be- 
tween the observed and expected Sl values for this model 
can be in part explained by the interplay of dissipational 
and dissipationless galaxy merging. Observations support 
the formation of more and metal-rich clusters during gas- 
rich mergers (e.g. Whitmore & Schweizer 1995; Whitmore 
et al. 1999; Goudfrooij et al. 2007). Semi-analytical models 
of cluster formation are able to reproduce the GC metallic- 
ity and age distributions assuming formation of GCs during 
galaxy merging (e.g. Muratov & Gnedin 2010). 

4.3 The Distribution of Globular Cluster 
Formation Efficiencies 

In the following we investigate the specific GC formation effi- 
ciency n as a function of galaxy luminosity and morphology. 
To derive the distribution of r) we choose Sm and Sl (upper 
panels in Fig. 6) as their values are least affected by internal 
systematics. Both parameters are independent of distance 
and are nearly insensitive to incompleteness of low-mass, 
faint GCs (see Sect. 3.2 and also Harris 1991, 2001). This is 
reflected in a smaller observed scatter in Sm and Sl values 
around the best-fit iso-77 relation compared to the Sn and T 
values. The advantage of the Sm value is that variations in 
stellar population parameters of the host galaxies and GC 
systems (i.e., their star-formation histories) are taken into 
account by the galaxy M/L (see Sect. 3.2). However, we 
point out that the derived Ty for low-mass dwarf irregular 
galaxies, which is based on galaxy colours, is likely to have 
an increased uncertainty due to their likely different SFHs 
relative to the sample of more massive galaxies from which 
the Bell et al. relations were derived (2MASS & SDSS galax- 
ies). We invert Equations (24) and (25) to derive rj values 
for our sample galaxies. The resulting rj distributions as a 
function of galaxy luminosity are shown in Figure 7. To test 
for trends of the GC formation efficiency as a function of 
galaxy luminosity and morphology we formed the running 
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average for rjL and t)m as a function of galaxy luminosity in 
galaxy luminosity bins of Mv = 1 mag, same as in Figure 4. 
The mean values are shown with solid curves and different 
symbols for the two galaxy morphology classes. The errors 
in the running average, not shown for clarity purposes, cor- 
respond to approximately d ex - 

The direct comparison between the 7]l and r/M distribu- 
tions in Figure 7 shows good agreement for early-type galax- 
ies which indicates a robust understanding of their mass-to- 
light ratios. However, for late-type galaxies the systemati- 
cally lower rj M with decreasing galaxy luminosity is caused 
by dwarfs without GC candidates. Those GC-free dwarf 
galaxies contribute mostly stellar light/mass but no GCs, 
i.e. the normalized GC formation efficiency is reduced at a 
given galaxy magnitude bin. 

Assuming consistent fj M and fj L distributions for dlrrs 
this translates into various (degenerate) systematics that are 
not included in the models, namely i) a growing mass com- 
ponent towards fainter galaxy luminosities, w) a steeper IMF 
in low-L dlrr galaxies, or/and in) a systematically higher 
internal extinction value in low-mass galaxies. The investi- 
gation of these systematics suggests a mass-dependent term 
a(Mh) m the scaling relation M h /L v <xM~ h 2/3+a(Mh) that 
is currently unaccounted for in the Dekel & Birnboim mod- 
els. The quantification of such a term goes beyond the scope 
of this paper and requires much more comprehensive GCS 
samples across the entire range of late-type galaxies from 
low-mass dwarf irregulars to massive spirals. 

The fits to the Sl — Mv and Sm — Mv distributions 
return (ij L ) = 5.5 x fCT 5 and (r) M ) = 6.7 x fCT 5 for the 
luminosity- and mass-normalized GC formation efficiencies, 
respectively. Those are shown with solid curves in the top 
panels in Figure 6. Because of the uncertainty in the derived 
galactic Ty affecting t)m we adopted the t]l = 5.5 x fCP 5 
for the bottom panels in Figure 6. The data can be approxi- 
mated by a wide range of specific GC formation efficiencies 
from 77 ~ 10~ 5 to ~ 5 x 10~ 4 . The varying scatter in rjL as 
a function of galaxy luminosity is at least partly due to the 
increasingly stochastic nature of star and star-cluster for- 
mation processes towards lower galaxy masses. This is sup- 
ported by the observational errors which are lower at higher 
galaxy luminosities (cf. Fig. 2) which relate to smaller stan- 
dard deviations of rjL (cf. top left panel of Fig. 7). At the 
faintest galaxy luminosities the jjl distribution is affected by 
the minimum-??!, relation (discussed in Sect. 4.f) which nar- 
rows its distribution. The large scatter in r/M at faint galaxy 
luminosities (relative to that of t)l) is primarily driven by 
the poor knowledge of the galaxy's MIL— ratios based on 
optical colours. 

The value of the mean GC formation efficiency, 77, that 
we obtain is consistent with previous estimates by Blakeslee 
(1999); McLaughlin (1999); Kravtsov & Gncdin (2005); 
Spitler & Forbes (2009). The value of r? 6 = 1.71 x 10~ 4 
obtained by Blakeslee (1999) is the GC formation efficien- 
cies with respect the baryon galaxy mass (stellar plus hot 
gas mass) derived for massive elliptical galaxies. Multiplying 
this value by fb = 0.17, the universal cosmological baryon 
fraction observed by WMAP (Hinshaw et al. 2009), leads 
to r)h = 2.9 x 10~ 5 . This value is identical to the 77— values 
obtained by us, Kravtsov & Gnedin (2005) and Spitler & 
Forbes (2009), rj h = 2 - 5 x 10~ 5 and 7.1 x 10~ 5 , respec- 
tively. 



To investigate how 77 relates to different galaxy mor- 
phologies we show histogram distributions in Figure 8 for the 
different morphological types: dlrrs, dSphs, dEs, spirals, and 
ellipticals. For those samples we derive the highest probabil- 
ity density distributions excluding the dwarfs with no blue 
or any GC detection (discussed in Fig. 7) . The number of 
such galaxies in the dlrrs, dSphs, dEs, Es and spirals sam- 
ples, respectively, is: 44 (51%), 4 (60%), 16 (18%), and 3 
(15%). For t]l, we find no evidence for a statistically sig- 
nificant offset in the average values (tjl) between late-type 
and early-type galaxies (see legend in Fig. 8) even though 
the scatter around the average value is larger for the fainter 
dwarf galaxies as mentioned above. We do however find a 
marginally significant offset between the (t)l) values of dlrrs 
and dEs. The fact that this offset is larger for (t)m) than for 
(t]l) suggests that it is mainly due to a systematic difference 
in stellar population properties between dlrrs and dEs. 



In general, an offset in rjL versus t)m between late- and 
early-type galaxies can be expected due to a difference in 
stellar population parameters of GCs and their host galax- 
ies. An increase in 7y by ~ 10%, which corresponds to a few 
Gyr age difference within the age range of our sample selec- 
tion (£gc > 4 Gyr) , introduces a shift in the rjL distribution 
towards larger values (Alogr/L = 0.03), which is insufficient 
to account for the observed difference in r\. However, lumi- 
nosity fading of the host galaxy of the order of AMv = 1.0 
to 1.5 mag (see Sect. 3.1) would cause log77L to shift by 0.24 
to 0.36 dex towards higher values, which would be enough 
to equalize the rjL distributions of late-type and early-type 
galaxies. Although some of our sample dlrrs may require 
such luminosity corrections (in particular those with very 
low 7]l), most of our sample dlrrs lack any evidence for sig- 
nificant present-day star formation rates which suggests that 
evolutionary fading corrections of their galaxy light might 
not be required. On the other hand, the mass-normalized 
GC formation efficiency rju should in principle correct for 
any stellar population age difference between late-and early- 
type galaxies. What we observe instead in Figure 8 is that 
the late-type galaxies in our sample (dlrrs and spirals) have 
lower average r\M than the early-type galaxies at the same 
luminosity, suggesting that the M / L— ratios of the late- type 
galaxies may be underestimated (as already discussed above 
in the context of Fig. 7). However, the number of spiral 
galaxies in our sample is only marginally significant to draw 
statistically meaningful conclusions. Therefore, a more com- 
prehensive study including near-infrared imaging or spectro- 
scopic observations (lacking for the majority of the dwarfs 
in our sample) should help in resolving this issue. 



Finally, we note that the minimum 77-limit described by 
Eq. (5) intercepts the iso-77 line for the maximum observed 
specific GC formation efficiency (Eq. 24) at Mv ~ —9 mag, 
corresponding to 3.37 x lO 5 !/© (cf. Fig. 6). In the context of 
our assumptions, this magnitude describes the lowest lumi- 
nosity halo that, on average, can form and host one old GC. 
This is consistent with the fact that the recently found ex- 
tremely faint Local Group dwarf galaxies do not possess GCs 
(e.g. Belokurov et al. 2006, 2007; Irwin et al. 2007; Walsh 
et al. 2007; Willman et al. 2005a,b; Zucker et al. 2006a,b). 
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Figure 7. Specific globular cluster formation efficiencies, 77, as a function of galaxy luminosity. Left: Luminosity normalized GC formation 
efficiency (t]l,) derived from the Sl relation (eqn. 24). Right: Mass normalized GC formation efficiency (vm) derived from the Sm relation 
(eqn. 24). The symbols connected with lines show the mean 77 values (including galaxies without GCCs) in bins of one magnitude. The 
remaining symbols are as in Figure 6. 



4.4 Implications for Hierarchical Galaxy 
Formation 

4-4-1 The Influence of various Physical Processes on the 
Variance of Globular Cluster Formation Efficiencies 

The general behavior of the GC system scaling relations im- 
plies that the star and star-cluster formation processes scale 
with the total host galaxy mass. However, it is not yet clear 
what is driving the large observed galaxy-to-galaxy scatter 
in the rj values for galaxies with otherwise similar mass. En- 
vironmental density is one of the suggested parameters to 
cause systematic trends in this regard (Peng et al. 2008). 
The majority of our current dlrr galaxy sample consists of 
field dlrrs, with a dearth of dlrrs situated in cluster environ- 
ments. Conversely, dE galaxies with GC data available in 
our literature compilation are all located in the Virgo clus- 
ter Miller & Lotz (2007); Peng et al. (2008). Hence we cannot 
yet make a statistically significant comparison between the 
GC formation efficiency of field and cluster dwarf galaxies. 
To robustly assess the impact of environmental density on 
the rj distribution will require the investigation of more dlrr 
galaxies in dense groups and galaxy clusters and more dE 
galaxies in loose environments. 

Another explanation for the scatter in 77 can be ob- 
tained from the different timing of star cluster vs. field star 
formation as a function of the star formation rate of the 
host galaxy. For the latter, Larsen & Richtler (2000) demon- 
strated that the specific [/-band luminosity in young massive 
star clusters in late-type galaxies strongly correlates with the 
total star formation rate and star formation rate surface den- 
sity of the host galaxy, SFR and Esfr, respectively. Thus, 
galaxies with the same mass but different SFR histories will 
have different GC formation efficiencies, since galaxies with 
higher peak SFRs will produce larger numbers of massive 
GCs that will survive a Hubble time. 

Furthermore, a delay between the peak formation 
epochs of star clusters and field stars that scales with SFR 
(in the sense that the star clusters form earlier on average 
than the field stellar population) will increase the values of 
the GC system scaling parameters, especially for more vio- 



lent starbursts which likely form more star clusters per unit 
gas mass. Any truncation of star-formation processes be- 
fore the formation of the field stellar population is complete 
(e.g., through gas stripping during infall into a dense galaxy 
cluster medium) would further increase the GC scaling pa- 
rameters. A hint that such effects might be at work was 
provided by Peng et al. (2008) who correlated the SFR and 
Ssfr with the mass fraction of GCs in their sample galax- 
ies to find that the cluster formation rate peaks ~ 500 Myr 
earlier than the total galaxy SFR. 

Another plausible mechanism to account for an appar- 
ent variation in 77 among galaxies is GC stripping and accre- 
tion (Cote et al. 1998, 2002; Hilker et al. 1999) that can cause 
a steeper slope of the Sl — Mv relation at brighter galaxy lu- 
minosities, where the accretion of high-77 dwarf galaxies may 
be an important mechanism during the early epochs of mas- 
sive galaxy assembly. Such dwarfs may have completed the 
process of star-cluster formation at high redshifts (2 >> 2) 
and had their field star-formation processes suppressed or 
stopped (due to shock heating) during their infall through 
the dense intracluster medium. 

4-4-2 Predictions from a Simplistic Accretion Model 

Some of the most luminous early-type galaxies in our sam- 
ple show excess 77 values (log rj « —3.5) that may indicate 
that their GCSs either experienced an early, highly-efficient 
episode of massive cluster formation followed by a shutoff 
of (field) star formation or that their GCSs were supple- 
mented with high-77 low-mass dwarf galaxies and their GC 
systems (or a combination of both effects). Such low-mass 
dwarfs would mainly contribute to the blue GC population 
and rather less to the field stellar component, and could well 
account for the observed high ratio of the metal-poor GCs to 
metal-poor field (halo) stars in massive galaxies (e.g., Harris 
& Harris 2002). Cote et al. (1998) and Hilker et al. (1999) 
suggested that faint (My « —12 mag) dwarfs accreted by 
cD galaxies can indeed provide the right number of GCs to 
account for their high Sm- We can make a simple estimate 
of the number and luminosity of satellite galaxies necessary 
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Figure 8. Specific GC formation efficiency, 77, for various galaxy 
types: dlrr, dSph, dE, elliptical and spiral galaxies (top to bot- 
tom panels). Left panels: luminosity normalized GC formation 
efficiency (vl) distributions. Those for dlrrs are not corrected 
for evolutionary fading of the integrated galaxy light. Right pan- 
els: mass normalized GC formation efficiency (r/jv/ ) distributions. 
Dashcd-dotted curves with Poisson errors are Epancchnikov- 
kernel probability density estimates. The maximum-likelihood r\ 
values for the different galaxy samples are shown in the upper left 
corner of each panel. The vertical solid line in all panels shows for 
comparison the 77— value for Es. Fading of AMy = 1.5 — lmag 
will shift the log n distribution of late-types toward higher values 
by 0.36 - 0.24 dex. 



to increase the Sl value of a progenitor galaxy by a certain 
value by simply using the definition of Sl (Eq. 3): 

-,2 iGCS.prg + A a cc X Z/GCS,acc 
Lv,pi S + A?acc x Lv, acc 



S L = lO 2 



(27) 



The above equation takes into account the increase in lumi- 
nosity of the progenitor galaxy and luminosity of its GCS as 
a consequence of the accretion. Assuming that the accretion 
process occurs in a dissipationless way (i.e. no new stars or 
GCs are formed) and that no GCs are destroyed, the num- 
ber of accreted galaxies can be expressed with the following 
relation: 



JVacc = 



A Si. 
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(28) 
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Figure 9. An example case for the number of satellite galaxies 
necessary to boost the specific luminosity of a host galaxy by 
ASl = 0-5 as a function of satellite luminosity, My acc . The 
relations are parametrized by the luminosity of the progenitor 
galaxy, My p , r The different line- types show three different cases 
of GC formation efficiencies of the accreted satellite galaxy (r? acc ). 
The peaks and "knees" in the curves are caused by the inflection 
point of the "U-shaped" GC system scaling relations (cf. Fig. 6). 



where Sl v _ prs , iv, P r E and Sl Viu . c , Lv,&cc are the progenitor 
and accreted galaxy specific luminosity and brightness, re- 
spectively. The increase in the GC system specific luminos- 
ity due to dissipationless accretion of satellite galaxies is 
ASl v = Sl,v — Sl Vpis - Sl,v is the specific luminosity 
of the remnant galaxy after the accretion processes ended. 
Relation 28 gives the prescription to compute the number 
of accreted satellite galaxies with a given luminosity and 
rja.ee (cf. top left panel of Figure 6) necessary to increase 
Sl by ASl- To illustrate relation (28), we show in Figure 9 
an example case for a progenitor galaxy with increased Sl 
value by ASl = 0.5 for three different 7; acc values of the 
satellite galaxies and progenitor luminosities in the range 
^v,pr g = —16 to —22 mag. The shape and slope of the 
curves is a direct reflection of the functional relation between 
the Sl value and galaxy luminosity observed in Figure 6 and 
described by eq. (24) (see also Eq. A28). The only purpose 
of this simple example is to show that in the formalism of 
the current model an accretion of low-luminosity dwarfs with 
high Sl values can increase the Sl of a more massive progen- 
itor galaxy. However, only a more detailed treatment with 
sampling an evolving galaxy luminosity function will yield 
realistic number estimates of accreted satellite galaxies dis- 
tributed over a mass spectrum, which is beyond the scope 
of this study. 

Thus, in this simplistic merging picture the accretion 
of a reasonable number of intermediate-luminosity, high-St 
dwarf galaxies offers a valid mechanism for boosting the GC 
system scaling parameters of cD galaxies. The observed cor- 
relation between the Sn of the cD galaxy and the X— ray 
flux of the host cluster Blakeslee (1999) supports this pic- 
ture. Accretion and galaxy merging in the cluster core heats 
up the gas leading to the high A— ray flux. Such accretion 
and hierarchical growth of cD galaxies is also evidenced 
by the steepening of the faint-end slope of the luminosity 
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function of galaxies in galaxy clusters with increasing red- 
shift (e.g., Ryan et al. 2007; Khochfar et al. 2007). In ad- 
dition, the fraction of low-mass peculiar/irregular galaxies 
increases with redshift (e.g., Renzini 2006; Conselice et al. 
2008). Moreover, observations of low-redshift Abell galaxy 
clusters show that galaxy clusters which host cD galax- 
ies have smaller dwarf-to-giant ratios than non-cD clusters 
(Lopez-Cruz et al. 1997) and a clear radial dependence of the 
faint-end slope of the galaxy luminosity function becoming 
steeper in the cluster outskirts (Barkhouse et al. 2007). It 
therefore seems that the contribution of low-mass dwarfs to 
the buildup of the GC system of the most massive galaxies 
is likely non-negligible, especially at high redshifts. In fact, 
stripping of GCs from present-day dwarf galaxies in the clus- 
ter environment is suggested by the lack of any detectable 
GC systems around few dEs closest to M 87 and M 49 in the 
Virgo galaxy cluster (Peng et al. 2008). In the context of the 
current cosmological ACDM paradigm of structure forma- 
tion (e.g., White & Rees 1978), it is of great importance to 
quantify the accreted mass fractions of GCs and field stars. 
It would be highly desirable to study the formation and en- 
richment timescales of the low and high-?; dwarf galaxies in 
more detail and compare the field stellar populations in such 
galaxies with those of their GC systems. 



5 CONCLUSIONS 

We investigate the much debated behavior of the observed 
GCS scaling parameters as a function of galaxy luminosity, 
such as the GC specific frequency (Sn), specific luminosity 
(Sl), specific mass (Sm), and specific number (T). Those 
are the integrated number, luminosity mass and specific 
number of all globular clusters in a galaxy normalized to 
the total galaxy luminosity and/or mass, respectively (see 
Sect. 3). We derive these quantities from HST/ 'ACS data 
of low-mass, faint (My > — 16 mag) dwarf galaxies, mainly 
late-type irregulars, located in the halo regions of nearby 
(< lOMpc) galaxy groups and in the field (Georgiev et al. 
2008, 2009a). In order to investigate the scaling relations 
of their GCSs as a function of galaxy luminosity (mass) we 
also compiled data from the literature for massive cluster el- 
lipticals (Peng et al. 2008) and spiral galaxies (Spitler et al. 
2008). To complement our data in the low-mass regime we 
included 69 cluster dEs (Miller & Lotz 2007) and 12 dlrrs 
from Seth et al. (2004) and Georgiev et al. (2006), as well 
as 24 late-type dwarfs in the local low-density environment 
from Sharina et al. (2005). Thus, we cover virtually the 
entire range in galaxy luminosity from My = — 11 to —23 
mag (1O 6 -1O 11 L ) where old GCs (i>4 Gyr) are reliably 
detected. Our main results can be summarized as follows: 

• Trends in GCS scaling relations hold irrespective of 
galaxy morphology. — The significantly increased number 
of galaxies in our analysis in comparison with earlier studies 
over the entire mass and galaxy morphology range allowed 
us to firmly corroborate the previous observational findings 
that the GCS scaling parameters vary as a function of 
galaxy luminosity (Fig. 6). We find that this relation holds 
irrespective of galaxy morphological type, suggesting a 
universal mode of GC formation. Galaxies show increasing 
GCS scaling parameters toward low and high-luminosity 



systems with a minimum at around Mv ~ —20.5 mag 
(L v « 2 x 1O 1O L ). 

• Sl values of early-type galaxies are twice that of late-types 
at a given galaxy luminosity. — The specific luminosity 
Sl is the most robust scaling parameter and shows least 
scatter due to its independence of distance and weak sen- 
sitivity to completeness corrections at the faint end slope 
of the GC luminosity function (Sect. 3.2.2). For late-type 
galaxies, spirals and irregulars the Sl value is on average 
two times smaller at a given galaxy luminosity than that 
for early- type systems (cf. Fig. 4). This difference can be 
partially accounted for by the passive evolutionary fading 
of the integrated galaxy light (AV > lmag), provided no 
more GCs are formed or destroyed at the same time (cf. 
arrow in Fig. 4). As cluster dSphs are on average fainter 
and exhibit higher scaling parameters than field dlrrs, 
the above analysis supports the idea that dlrr galaxies 
could be the progenitors of dSphs and fainter dEs (e.g. 
Miller et al. 1998), but see also Grebel et al. (2003). 
However, our sample contains mostly dlrrs in field and 
group environments. Their GC scaling parameters might 
differ in cluster environments due to varying GC formation 
efficiencies and/or cluster stripping. However, such cluster 
dlrrs have not been sampled well enough yet. 

• The "U-shaped" behavior of the GC scaling parameters is 
remarkably well described by a two component model at a 
transitional halo mass critical for the thermal properties of 
the inflowing gas. — In order to explain the trends in the 
GCS scaling parameters we have assumed that the total 
mass in GCs is proportional to the total halo mass of the 
host galaxy (TMgcs = r/Mh), i.e. the formation of GCs 
scales with the total galaxy mass. The coefficient, rj, is the 
specific GC formation efficiency parameter, which measures 
the present-day mass/luminosity of the GCS drawn from 
the initial star cluster population mass/luminosity function 
shaped by various mechanisms over a Hubble time (cluster 
relaxation, tidal and dynamical friction, dissolution). We 
have invoked theoretical models of Dekel & Birnboim 
(2006) which predict the dependence of the galaxy M/L 
as a function of galaxy mass determined by the thermal 
properties of the inflowing gas for two different halo mass 
regimes in which star formation processes are regulated 
by supernova feedback and virial shocks below and above 
the critical stellar mass of 3 x 10 10 X Q . For low- mass 
halos these models predict that the mass-to-light ratio 
varies with host galaxy mass as Mh/L oc M^ 2 ^ (see also 
Dekel & Silk 1986). To draw quantitative conclusions, the 
absolute normalization of this relation is very important 
which we derive from recent observations of the mass in 
the inner 0.3 kpc (.Mo. 3 = 10 7 yV/fo) of nearby low-mass 
dwarf galaxies (Strigari et al. 2008) using its relation to the 
total halo mass from the models of Bullock et al. (2001). 
Above the critical mass, at large halo masses the models 
predict Mh/L oc M\[ 2 which we calibrate with dynamical 
mass measurements based on galactic group dynamics of 
massive early- type galaxies (Eke et al. 2006). This results 
in an analytical model which describes remarkably well the 
observed trends in the GC scaling relations as a function 
of galaxy luminosity (Sect 4.2.2). In the Dekel & Birnboim 
(2006) model, compact and dense molecular clouds (the 
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sites of GC formation) are efficiently shielded against shock 
heating and ionizing stellar feedback above and below the 
critical galaxy mass, respectively. This implies that in this 
picture, the shape of the GC scaling relations can be de- 
scribed as a universal GC formation efficiency for the entire 
galaxy mass range but with an evolving field star-formation 
efficiency. The latter is equally low for low- and high-mass 
systems (compared to the formation of massive/dense 
star clusters) due to being efficiently suppressed by the 
stellar and SNe ionizing feedback and shock heating of the 
inflowing gas below and above the critical galaxy mass 
(~fO lo ' 5 X ), respectively. This causes an evolution of the 
mass in clusters to mass in field stars ratio as a function 
of galaxy mass. The evolution of this ratio is at the roots 
of the shape of the GC scaling relations. The best fit of 
this model to the observed distributions yielded a value of 
the observed GC formation efficiency parameter (77) for the 
luminosity and mass-normalized GC formation efficiency of 
{r) L ) ^ 5.5 x 1(T B and (n M ) ^ 6.47 x 10" 6 (see Sect. 4.3). 
The 77 distributions are very similar for all galaxy types, if 
passive evolutionary fading of the late-type galaxy sample 
of ~ f mag is applied (cf. Fig. 8). 

• The differences between model predictions and observa- 
tions can be attributed to a mixture of effects of varying 
GC formation efficiencies, galaxy merging histories, and a 
variation in cluster destruction mechanisms as a function 
of galaxy mass lacking in the current model. — The most 
massive galaxies, whose Sm and Sl values increase more 
rapidly than expected from the theoretical predictions 
for a fixed 77 can be understood as systems which have 
either undergone an early episode of extremely efficient 
star-cluster formation or less efficient formation of field 
stars or which have preferentially accreted high-77 dwarf 
galaxies. In a simple merging picture of satellite galaxies 
of a given luminosity we showed in Sect. 4.4 that it is 
possible to boost the Sl values of galaxies by accretion 
of intermediate luminosity, high-77 dwarfs, thus offering an 
efficient mechanism for explaining the high 77 values of cD 
galaxies. In addition, the observed spread in 77 may be 
caused by the stochastic nature of star and star-cluster 
formation at low galaxy mass or by differences in SFH and 
SFR intensity for galaxies at the same luminosity. 

To understand whether the difference in the GCS 
scaling-parameter (Sl or Sm) distributions between late- 
and early-type galaxies is statistically significant and, if so, 
to understand the nature of this difference (e.g. age, chemical 
enrichment, environment, etc.) it is crucial to sample with 
more observations (e.g. near-UV/IR photometry or spec- 
troscopy) the GCSs of spiral galaxies at intermediate to high 
masses, late-type dwarfs in dense cluster environments, both 
of which are significantly under represented in the current 
study. 
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APPENDIX A: SCALING RELATIONS OF GC 
SYSTEMS AS A FUNCTION OF GALAXY 
MASS (DARK AND BARYON) AND 
LUMINOSITY 

Al GCSs scaling relations in the "SNe feedback" 
galaxy mass regime 

Using Equations (9) through (11) we express Equations (2), 
(3), and (4) as a function of galaxy luminosity, total halo 
and baryon mass for M h <3 x 1O 1O M . 



= 2 + log(^-) +logM b 
log S M = 2 + log [^^j + log L v 

= 2 + lo s(^) + log Mb 



log S N = 6 + O.4M y ,0 + log 



771TO 



0.4 log L v (Al) logT = 9 + log 



TyJTlTO 



+ log L V 



= 6 + O.4Mv, + log (^L) - § lo gMh (A2) = 9 + log + \ ^ Mh > 



6 + OAMvq + log I -5_ ) - 0.4 log M h (A3) 



log S L = 2 + log^^-^ -0.4 log L v 
= 2 + log(^)-|log^ 
= 2 + log P" 1 Iy j -0.41ogM b 

logSW = 2 + \og( T ^) -0.4 log Lv 



= 2 + log 



Tv 

0.6 



3 logX h 



S + logp^ ) -°- 41 °S- Mb 



v 

0.6 



logf ^9 + logf- 7 ^) 



= 9 + log 



TyriiTO 
0.6 



0.4 log Lv 
2 



9 + 1OS (^^)- - 41Og ^ 



(A4) 
(A5) 
(A6) 
(A7) 
(A8) 
(A9) 
(A10) 
(All) 
(A12) 



A2 GCSs scaling relations in the " Virial-shock" 
galaxy mass regime 

Using Equations (9) through (12) we express Equations (2), 
(3), and (4) as a function of galaxy luminosity, total halo 
and baryon mass for Mh > 3 x lO 10 ^©. 

log S N = 6 + 0.4Mv © + log ( \ + 1.2 log L v (A13) 

\m T o / 

= 6 + O.4Mv, +logf^^ N ) +hogM h (A14) 



= 6 + O.4Mv, +log 



+ 1.21ogM b (A15) 



log SL = 2 + log (^^j + log Lv 

= 2 + logf ?7K2N ) +^logM h 
V 7v J 2 



(A16) 
(A17) 



9- loo ( -5-^- ] H- Ior.M,,. 
m T o T v 



(A18) 
(A19) 
(A20) 
(A21) 
(A22) 
(A23) 
(A24) 



A3 Continuous GCSs scaling relations from the 
"Virial-shock" to the "SNe feedback" galaxy 
mass regime 



S N = 10 6+0 - 4M ^ -5_ ( K ?- 6 Lv°' 4 + 4L V 2 ) 
m T o 



(A25) 



1O 6+O.4M V , _V_ ( KlM -W + ^°- 8 M 3 h /5 ) (A26) 
m,TO ^ ' 

= 10 6 + ' 4M ^ ^- ( ( Kl Ty)°- 6 M b - 4 + ^M b 2 ) (A27) 



c , n 2 1) ( 0.6 r -0.4 . 2 r \ 

Sl = 10 — Ki L v +k 2 L v ) 

7v v ' 

= io 2 ^L( KlA1 - 2/3 + ^ /2 ) 

7v V / 



= io 2 ^ ur?VM^+^M b 

7v V 1 v 

c , n 2 1 / 0.6 r -0.4 . 2 r \ 

Sm = 10 — (Ki L v + n 2 L v ) 



= 1{ ?JL( Kl M?'* + 



T= 10 



9 »7 



/ 0.6 r -0.4 . 2 r \ 
(Kl L v +K 2 Lv) 



10 9 ^V( Kl M- 2/3 + K2 A< /2 ) 



10 



9 f] 



K2_ 

m T o \ \Tv 



( f ; ) m, 



(A28) 
(A29) 

(A30) 
(A31) 
(A32) 

(A33) 

(A34) 
(A35) 

(A36) 



